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Abstract. We explore the feasibility of estimating primary cosmic ray 
composition at high energies from the study of two parameters of Ex¬ 
tensive Air Showers (EAS) at ground and underground level with Monte 
Carlo simulations using the new EPOS and QGSJETII hadronic models 
tuned with LHC data. Namely, the slope and density at a given distance 
of the muon lateral distribution function are analysed in this work. The 
power to discriminate primary masses is quantified in terms of merit fac¬ 
tor for each parameter. The analysis considers three different primary 
particles (proton, iron and gamma), four different zenith angles (0°, 15°, 
30° and 45°) and primary energies of 10 17 ' 25 eV, IQ 17 - 50 eV and 10 17 75 eV. 


1. Introduction 

The particle lateral distribution of EAS is the key quantity for cosmic ray ground 
observations at energies greater than 10 15 eV, from which most observables are 
derived. An EAS is initiated by a high energy cosmic ray interacting in the top 
of the atmosphere and creating a multitude of secondary particles, which arrive 
at ground nearly at the same time. Secondary particles are distributed over a 
large area perpendicular to the direction of the cosmic ray primary. The disc 
of secondary particles may extend over several hundred meters from the shower 
axis, reaching its maximum density in the center of the disc, which is called 
the shower core. The density distribution of particles within the shower disc 
can be used to derive information on the primary particle. EAS measurements 
at ground level are carried out using arrays of individual detectors, which take 
samples of the shower disc at several distances from the shower core [1|. 

It is known that the muonic component of EAS carries information about 
the identity of the primary particle [2]. In this work we study two parameters 
sensitive to primary masses: the slope of the Muon Lateral Distribution Function 
(MLDF) and the density of muons at a certain distance from the core. The 
underlying idea is that showers originated by heavy nuclei produce more muons 
than lighter ones (therefore the absolute muon density at fixed distance will be 
higher for heavy nuclei) but the MLDF will be less steep. For the presented 
analysis the response of any detector is not taken into account. 
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2. Monte Carlo simulations 


For this work we generated a library of EAS using AIRES 2.8.4a [3] which make 
use of QGSJET-II-03 [3| as hadronic model and CORSIKA 7.3700 [5] which 
make use of QGSJET-II-04 @] and EPOS-LHC [ 6 ] as hadronic models. In both 
programs, we set a statistical thinning of 10 -6 . For each hadronic model we 
consider three types of primaries (proton, iron and gamma), four zenith angles 
(0°, 15°, 30° and 45°) and three energies (10 1725 eV, 10 1 '" 50 eV and 10 17 ’ 75 eV). 
For each energy, zenith angle and primary type a total of 120 showers were 
produced, considering an uniform azimuthal distribution between 0° and 360°. 
During the EAS simulations only muons with energies above 55 MeV have been 
taken into account to assess the MLDF at different depths between 0 and 2.5 m 
using the same set of showers. 


3. /3 and p M (500) parameters 

Once the EAS is simulated, the muon lateral density is fitted event-by-event with 
a KASCADE-Grande like MLDF [7] : 


p M (r) = Np 




(1) 


where r is core distance in the shower plane. The values of a, 7 and r 0 are fixed 
in 0.75, 3 and 320 m, respectively, and and j3 are the free parameters. 

The MLDF is simulated at ground with AIRES and CORSIKA and afterwards 
propagated underground. In the last case, two depths have been used: 1.3 m 
and 2.5 m. The layer of soil is used as a shielding against the electromagnetic 
component of the EAS, allowing only the muons of energy greater than 0.52 and 
1 GeV to arrive to the desired depths respectively. For each event simulated, /3 
and p ti (r = 500) are obtained from Eq. (jT]) . For the energies considered, the 
distance of 500 m is close to the distance where the fluctuations over the MLDF 
are minimized if the any detector response is not taken into account. However, 
Fig. [I] (Right) show a wide range in distance with similar fluctuations values. 


4. Muons propagation through the soil 

Muons lose a fraction of their energy when they propagate through the soil, 
mainly due to ionization. Therefore, we assume as a first-order approximation 
that the energy loss is proportional to the muon track length and constant with 
energy. Then the energy of a muon that traveled a distance x through the soil 
is: E /J/ (x) = Efj, 0 — apx, where 0 is the initial energy of the muon, p = 1.8 x 10 6 
g m -3 is the density of a soil with standard rock and a = 1.808 x 10~ 7 GeV m 2 
g -1 is the fractional energy loss per grammage depth. As an example, in Fig. 
|T]we show a fit of the muon lateral density for one vertical proton with energy 
of 10 17 ' 75 eV and QGSJET-II-03 hadronic model. In this case, the slope values 
were /3 = 2.44 ± 0.38, /3 = 2.67 ± 0.43 and /3 = 2.93 ± 0.49, and the density values 
were p M (500) = 0.81T0.12 nr” 2 , p M (500) = 0.62T0.09 m“ 2 , and p M (500) = 0.48T0.07 
m~ 2 , for ground, 1.3 m and 2.5 m underground, respectively. 

In order to validate the first-order approximation, we use the GEANT4 package 
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Figure 1. Left: MLDF of one vertical proton of 10 17 ' 7 ' 5 eV at different 
soil levels. The muons propagation was performed using the first-order 
approximation. Right: The mean spread in pn(r) for 120 showers initi¬ 
ated by proton with energy of 10 1 ' 50 eV and zenith angle of 30°. 

to perform a complete simulation of particles arriving to different depths. The 
results are equivalent because the ionization is the main process of energy loss. 


5. Data Analysis 

In order to quantify the discrimination power of (3 and p^(500) parameters, we 
constructed the distributions of these two parameters for each primary type, 
energy and zenith angle, at ground and underground levels, respectively. As 
an example, in Fig. [2] the distributions of (3 and p^(500) are shown for E = 
10 1 '' 75 eV and zenith angle 0°. The power of discrimination is measured cal¬ 
culating the Merit Factor (MF) between the distributions defined as: MF = 
l<A> - (B)\/^/a 2 A +a%, where (A) and ( B ) are the mean values of the A and B 
distributions, and a a and (Jb their respective standard deviations (because no 
detector is considered in this work, these fluctuations correspond only to shower- 
to-shower fluctuations). MFs lower than 1 indicate very poor discrimination 
power (roughly because mean values are less than 1 o apart from each other). 



Figure 2. Left: (3 parameter distributions corresponding to 120 verti¬ 
cal EAS at ground and 2.5 m underground levels, with energy of 10 17 ' 5 
eV and QGSJET-II-03 hadronic model. Right: The same for p^{ 500) 
parameter distributions. 
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6. Merit factor and primary mass discrimination 


The MFs dependence with energy, hadronic model and zenith angle was the main 
goal of this work. In Fig. [3] we show the results for case of 2.5 m underground. 
From both plots it can be seen that the discrimination power of (3 is not very 
good. On the other hand, p^(500) shows a MF greater than 2 and is practically 
independent of the energy and zenith angle. Therefore this parameter has the 
potential to discriminate the mass of primaries cosmic rays. A similar behavior 
was seen at 1.3 m underground level. 


Gamma and Proton 


Proton and Iron 




Figure 3. Merit factor of the ^(500) and (3 parameter at 2.5 m un¬ 
derground level as function of primary energy and zenith angle. Left: 
For Gamma and Proton. Right: For Proton and Iron. 

7. Conclusions 


In this work we performed a study of the merit factor of [3 and ^(500) pa¬ 
rameters obtained from MLDFs. We used proton, iron and gammas as primary 
cosmic ray, three hadronic models and different ground levels. The muons prop¬ 
agation through of soil was performed using an approximation which only took 
into account the continual energy loss. For energies and zenith angles studied, 
the f3 parameter has a merit factor < 1, and therefore it is not a good mass 
discriminator, while the ^(500) parameter has a merit factor > 2 (if Poissonian 
fluctuations are not taken into account). In each case, the merit factor does not 
show a significant difference with the three hadronic models considered. 
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